Fault-tolerant holonomic quantum computation 
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We explain how to combine holonomic quantum computation (HQC) with fault tolerant quantum 
error correction. This establishes the scalability of HQC, putting it on equal footing with other mod- 
els of computation, while retaining the inherent robustness the method derives from its geometric 
nature. 
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Introduction. — Holonomic quantum computation 
(HQC) 1] is an all-geometric method of computation 
that makes use of non-Abelian generalizations of the 
Berry phase 0. It is a fundamental alternative to 
other established models of quantum computation (QC) 
which in some sense lies between the circuit and 
adiabatic models, in that it replaces the dynamical 
gates of the circuit model by adiabatic holonomies. In 
HQC states are encoded in the degenerate eigenspace of 
a Hamiltonian and gates are realized by adiabatically 
varying the Hamiltonian along suitable paths in pa- 
rameter space, giving rise to geometric transformations. 
HQC has attracted significant interest both because of 
its deep connections to gauge theory Q and its potential 
practical advantages. It has been shown that due to its 
geometric nature, the method is resilient to various types 
of errors in the control parameters driving the evolution 
Q and thus could provide a certain built-in robustness 
at the hardware level. On the other hand, geometric 
phases are susceptible to decoherence, a conclusion 
which also affects HQC 0. 

No method of computation is scalable without the abil- 
ity to implement it fault-tolerantly, and HQC is no ex- 
ception. However, unlike the circuit and one-way mod- 
els of quantum computation Q for which fault tolerance 
proofs exist under reasonable noise models @, @, [§], a 
demonstration that HQC can be made fault tolerant un- 
der similar assumptions has been lacking. Here we rem- 
edy this situation, establishing the in-principle scalability 
of HQC. 

Fault-tolerant techniques guarantee that, under appro- 
priate assumptions such as errors that are sufficiently 
uncorrelated and improbable, an arbitrarily large com- 
putational task can be implemented with relatively mod- 
est resource overhead that preserves the computational 
complexity class. This result, known as the threshold 
theorem, is based on the use of quantum error correcting 
codes (QECCs) — a general software solution to the 
problems of noise and decoherence in quantum comput- 
ers. In Ref. [llj], a fi rs t step was taken to protect HQC 
against decoherence, by combining it with the method of 
decoherence- free subspaces (DFSs) leading to pas- 
sive protection against certain types of correlated errors. 



However, this is not enough for fault tolerance, since 
other error types could accumulate detrimentally unless 
corrected. Therefore, scalability of HQC requires going 
beyond the scheme of Ref. [HJ, e.g., by combining the 
holonomic approach with active error correction. Here 
we present a scheme for fault-tolerant HQC using sta- 
bilizer QECCs We begin by briefly reviewing the 
basics of HQC and the main ingredients of fault-tolerant 
QC. We then show how these ingredients can be realized 
using holonomic transformations on QECCs, and prove 
that the construction is fault-tolerant. Finally, we dis- 
cuss the properties of the scheme and analyze its effect 
on the accuracy threshold. 

Holonomic quantum computation. — Let {^a} be 
an isodegenerate family of Hamiltonians on an N- 
dimensional Hilbert space, which is continuously pa- 
rameterized by a point A in a differentiable control- 
parameter manifold M, i.e., H\ = J2n=i £ n(X)H n (X), 
where {e n (A)}^ =1 are the R different d„-fold degen- 
erate eigenvalues of H\, (%2 n=1 d n — N), and n n (A) 
are the projectors on the corresponding eigenspaces. If 
the parameter A is changed adiabatically (i.e., suffi- 
ciently slowly to prevent transitions between different 
eigenspaces), the unitary evolution that results from the 
action of the Hamiltonian H(t) := H x m is 

U(t) = Texp(-i [ drH(r))=®^ 1 e i ^U^ n (t), (1) 
Jo 

where oj n (t) — — L rfT£„(A(r)) is a dynamical phase, and 



,A(t) 

Ui l (t)=Vexp{ Y^A^dXn 



(2) 



Here T and V denote time- and path-ordering. The 
adiabatic connection is ^ M An^dX 11 , where A M are lo- 
cal coordinates on At (1 < /i < dinxM), A n ^ 
has matrix elements [1] (A nifi ) a p — (not; A|^jjr|n/3; A), 

and {\na; X)}'^L 1 is an orthonormal basis of the n th 
eigenspace of the Hamiltonian at the point A. 

When the path A(i) forms a loop j(t), j(0) = 
7(T) = A , the unitary matrix U% = U\ n {T) = 
Vey^(§ ^2 An^dX^) is called the holonomy associated 
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with the loop. When the n energy level is non- 
degenerate (d n = 1), the corresponding holonomy re- 
duces to the Berry phase 0] . The space of all loops based 
on A is L Xo (M) = { 7 : [0,T] M\j(0) = 7 (T) = A }. 
The set Hol(A„) = {U^Yi € £a (A4)} is a subgroup of 
U(d n ) called the holonomy group. If the dimension of M 
is sufficiently large, non-Abelian holonomies can be used 
to implement universal quantum computation over states 
encoded in one of the degenerate eigenspaces of H(t) 0. 

Fault-tolerant operations. — We are concerned with 
standard [l3T ] and operator stabilizer QECCs for 

the correction of single-qubit errors, and the techniques 
for fault-tolerant computation 0, H| on such codes. A 
quantum information processing scheme is called fault- 
tolerant if a single error occurring during the implemen- 
tation of any operation introduces at most one error per 
block of the code (a block is a set of qubits encoding one 
logical qubit). It is known 7] that fault-tolerant infor- 
mation processing is possible on any stabilizer code and 
can be realized almost entirely in terms of transversal 
operations — these are operations for which each qubit in 
a block interacts only with the corresponding qubit from 
another block or from a special ancillary state. In addi- 
tion, it is required that we are able to prepare and verify 
a special ancillary state. Since single-qubit unitaries to- 
gether with the "controlled-not" (CNOT) gate form a 
universal set of gates, fault-tolerant computation can be 
realized entirely in terms of single-qubit operations and 
CNOT operations between qubits from different blocks, 
assuming that the special state can be prepared reliably. 
Hence, our goal is to construct holonomic realizations of 
these operations, as well as of the operations needed for 
the preparation and use of the ancillary state. 

The manner in which we combine HQC and QECCs 
is by embedding the entire stabilizer code space or code 
subsystem into the ground space of a two-level degener- 
ate Hamiltonian H\ or in a subsystem which is invariant 
under the action of H\. The Hamiltonian is an element 
of either the stabilizer or the gauge group of the code at 
the initial moment and remains an element of the trans- 
formed stabilizer or gauge group at every moment dur- 
ing the computation. We perform computation in both 
the ground and excited eigenspaces of H\. Note that 
in this regard we depart from the original HQC method 
[l|, where computation is performed entirely within a 
single eigenspace. Even though the geometric approach 
requires the use of degenerate Hamiltonians which un- 
avoidably couple qubits within the same block, we show 
that propagation of single-qubit errors can be avoided. 

The scheme. — Consider a stabilizer code for the correc- 
tion of arbitrary single-qubit errors. We first show how 
to implement encoded Clifford gates on such a code (Clif- 
ford gates are those that preserve the Pauli group — the 
group of tensor products of Pauli operators). It is known 
[7| that these gates can be realized using transversal Clif- 
ford operations. For simplicity we restrict our attention 



to implementing transversal operations on the first qubit 
in each block; the operations on the remaining qubits 
can be obtained analogously, and used to complete the 
encoded Clifford gates on or between code blocks. As 
a starting Hamiltonian for implementing a single-qubit 
operation, we choose an operator that is an element of 
the stabilizer, or the gauge group (for the case of sub- 
system codes [IH), and acts non-trivially on that qubit. 
Without loss of generality we can write the initial Hamil- 
tonian as H(0) = —Z <g> G, where X, Y and Z are the 
Pauli matrices, and G is a tensor product of Pauli ma- 
trices on the rest of the qubits in the block and possibly 
on qubits from other blocks if we are in the middle of an 
entangling operation. (We can assume that G spreads 
over at most 4 blocks, since this is sufficient for imple- 
menting transversally any encoded Clifford operation on 
stabilizer codes [7j.) Henceforth a hat denotes operators 
on all qubits and tilde on all qubits excluding the first 
one. It is not hard to show [l6| that if the Hamiltonian is 
varied adiabatically so that only the first factor changes, 

H(t) = -H(t)®G, (3) 

where Tr{iJ(t)} = 0, then the geometric transformation 
resulting in each of the eigenspaces is equal to a local 
unitary operation on the first qubit: 

U(t) =U(t)® 7; U(t) = u Ao (t) e U Al (t) , (4) 

where U A Jt) =e& dr ^"WI^I^W>|^ n (t))(0 n (Q)|, with 
\4>n(t)} being the ground (n — 0) and excited (n = 1) 
states of H{t). 

Equation ((4]) is the basis of our construction. If H(0) 
is minus an element of the stabilizer, then the code space 
belongs to its ground space. Assuming that the encoded 
state has not undergone an error, by varying the factor 
H(t) adiabatically we can effectively generate any single- 
qubit unitary transformation on the state (we show how 
to do this below). If the initial Hamiltonian is an oper- 
ator in the gauge group for the case of subsystem codes, 
the non-erroneous state of the system can be a superposi- 
tion of ground and excited states. According to Eq. (J4|, 
each of the ground and excited components will undergo 
the same single-qubit unitary transformation U(t), but in 
addition, a relative phase of dynamical origin will accu- 
mulate between the two. This relative phase is equivalent 
to an operation on the transformed gauge subsystem, and 
therefore does not affect the encoded state. 

Single-qubit operations. — We now show how the 
method above can be used to gene rate a set of standard 
single-qubit gates. It turns out [16| that any Hamiltonian 
of the type © where 

H(t) = f(t)Z + g{t)V e ±ZV ] e± = f(t)Z + g(t)H g± (5) 

with f(T) = g(0) = 0, and f(0),g(T) > 0, gives rise in 
the adiabatic limit to the geometric single-qubit unitary 
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transformation 



V e ± = 



1 Te" 
1 



(6) 



where 8 is a real parameter. Let us define the eigen- 
states of H{t) at time T as \<p n (T)) = V 9 ±\n), n e {0, 1}. 
We can then write UA n (T) — e 1 "" Ve±\n)(n\, where a n 
are geometric phases, and we can show that e %a ° — e tai 
[ill]. Therefore, up to a global phase, Eq. (f4|) yields 
U(T) = Vg±. Using the last result and the identity 
±(cos9X + sin^F) = Hg±, one can verify the follow- 
ing constructions: an adiabatic interpolation along the 
path —Z® G — > —X ® G yields the operation Vq+ = RZ 
where R is the Hadamard gate; an adiabatic interpolation 
along the path-Z<g)G^ -(4=X + iY) (g) G — > Z®G 

yields %/4_^7r/4+ = V^^ 4+ , which up to an overall phase 
is equal to XS where S is the Phase gate. RZ and XS 
can generate all single-qubit operations in the Clifford 
group. Note that single-qubit gates outside the Clifford 
group can also be implemented in a similar manner [l6| . 
We point out that the change of the Hamiltonian along 
the edges of the paths must be sufficiently slow so that 
the adiabatic approximation 17[ is satisfied within the 
desired precision; we return to the validity of the adia- 
batic approximation below. 

Completing the gate set. — To complete the set of gates 
needed for encoded Clifford operations, we only have to 
show how to implement a transversal CNOT gate with 
the first qubit being the target. At the moments between 
the basic operations building up an encoded Clifford gate, 
we can always find G which acts trivially on the control 
qubit. Then the CNOT gate can be applied by first ap- 
plying the inverse of the Phase gate on the control, as 
well as the transformation —Z®G — ► —Y®G on the 
target, followed by the transformation — I c <g) Y <g> G — > 
— Z c ® Z® G where the superscript c denotes the control. 

Encoded operations outside of the Clifford group re- 
quire the above transformations plus the ability to mea- 
sure a particular encoded Clifford operator Q • The latter 
involves applying the operator conditioned on the qubits 
in a "cat" state (|0...0) + |l...l))/\/2. If we were im- 
plementing solely the Clifford operator, at any stage be- 
tween the elementary gates we would use a stabilizer or 
gauge group element of the form G — G\ ® Gj, where 
G\ is a tensor product of Pauli matrices acting on the 
first qubits from the blocks, and G\ is an operator on 
the rest of the qubits. Applying a basic transversal op- 
eration O conditioned on the first qubit in a cat state 
transforms this operator as I c <g> G\ ® Gj — > |0)(0| c (g> 
Gi®Gi + |l)(l| c ®OGiO t <£>Gi, where the superscript c 
denotes the control qubit from the cat state. The corre- 
sponding gate can be implemented via the Hamiltonian 

H (t) = -|0)(0| c ®Gi®Gi -a(i)|l)(l| c ®ff (i)®Gi, 
where Ho (t) <£> Gj is the Hamiltonian that we would use 
for the implementation of the operation O, and a{t) is 



chosen so that the operator a(i)|l)(l| c ® Ho it) ® Gj has 
the same instantaneous spectrum as |0)(0| c ® G\ ® G\. 
(Any possible relative geometric phase between |0) c and 
|l) c can be corrected by a suitable gate on the control 

qubit.) Starting from Ho(T), we can implement another 
conditional transversal operation in a similar fashion, etc. 

Fault tolerance of the scheme. — So far we have shown 
how to generate any transversal operation on an en- 
coded state holonomically, assuming that the state is 
non-erroneous. But what if an error occurs on one of 
the qubits? At any time t, we can distinguish two types 
of errors: those that result in transitions between the 
ground and excited spaces of H(t), and those that result 
in transformations inside the cigcnspaces. Due to the dis- 
cretization of errors in quantum error correction (QEC), 
it suffices to prove correctability for each type separately. 
The key property of our construction is that the geomet- 
ric transformation is the same in each of the eigenspaces, 
and it is transversal. Because of this, if we are applying 
a unitary on the first qubit, an error on that qubit will 
remain localized regardless of whether it causes an excita- 
tion or not. If the error occurs on one of the other qubits, 
at the end of the transformation the result would be the 
desired single-qubit unitary gate plus the error on the 
other qubit, which is correctable. It is remarkable that 
even though the Hamiltonian couples qubits within the 
same block, single-qubit errors do not propagate. This 
is because the coupling between the qubits amounts to 
a change in the relative phase between the ground and 
excited spaces, but the latter is irrelevant since it is ei- 
ther equivalent to a gauge transformation, or when we 
apply a correcting operation we project on one of the 
eigenspaces. In the case of the CNOT gate, an error can 
propagate between the control and the target qubits, but 
it never results in two errors within the same codeword. 

In addition to transversal operations, a complete fault- 
tolerant scheme requires the ability to prepare, verify 
and use special ancillary states, e.g., Shor's cat state [8j. 
Since the cat state is known and its construction is non- 
fault-tolerant, this can always be done using our holo- 
nomic approach by treating each initially prepared qubit 
as a simple code (with G being trivial) , and updating the 
stabilizer of the code via the applied geometric transfor- 
mation as the operation progresses. The only difference 
is in the measurement of the parity of the state, which re- 
quires the ability to apply successively CNOT operations 
from two different qubits in the cat state to one and the 
same ancillary qubit initially prepared in the state |0). 
After this operation, the ancilla qubit is measured in the 
{|0), |1)} basis so the relative phase between these two 
states is irrelevant. We can regard the qubit in state |0) 
as a simple code with stabilizer (Z), and we can apply 
the first CNOT as described before. Even though after 
this operation the state of the target qubit is unknown, 
the second CNOT can be applied via the same interac- 
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tion, since the transformation in each eigenspace of the 
Hamiltonian is the same. 

Discussion. — Since the method we presented conforms 
completely to a given fault-tolerant scheme, it does not 
affect the error threshold per operation for that scheme. 
Some of its features, however, affect the threshold for 
environment noise. First, observe that when applying 
the Hamiltonian ([3]), we cannot at the same time ap- 
ply operations on the other qubits on which the factor 
G acts non-trivially. This could be a problem for non- 
concatenation schemes, but with concatenated codes it 
only affects the lowest level of concatenation — some op- 
erations at the lowest level that would otherwise be im- 
plemented simultaneously might have to be implemented 
serially. This has the effect of slowing down the circuit 
by a small constant factor. For example, it turns out that 
for the 9-qubit Bacon-Shor (BS) subsystem code 14| this 
slowdown is by a factor of 1.5 [161 ]. 

A more significant slowdown results from the fact that 
the evolution is adiabatic (however, as argued in Ref. 
[Til ], an adiabatic implementation may be unavoidable 
in the case of Markovian noise, since fast gates are in- 
compatible with the Markovian limit). In order to ob- 
tain a rough estimate of the slowdown due specifically 
to the adiabatic requirement, we compare the time Th 
needed for the implementation of a holonomic gate with 
precision 1 — 6 to the time Td needed for a dynamical 
realization of the same gate with the same strength of 
the Hamiltonian. We consider a realization of the X 
gate via the Hamiltonian H(t) = V x {y{t))ZV}r(y{t)), 
V x (y(t)) = exp(iy(t)jfcX), where y(0) = 0, y{T h ) = T h . 
The energy gap of this Hamiltonian is constant. The op- 
timal dynamical implementation of the same gate is via 
the Hamiltonian —X for time Td = ^. It is known that 
if H(t) is smooth (but non-analytic) and its derivatives 
vanish at t = and t = Th, the adiabatic error decreases 
super-polynomially with Th [13] • To achieve this, we 
choose y(t) = ± f Q dt'e' 1 / siD( ^' / Th \ a = y(T h ). For this 
interpolation, when Th/Td ~ 17, the error 5 is already 
~ 1CD 6 , which is below the threshold values obtained, 
e.g., for the BS codes (~ 1CT 4 ) 0. This slowdown 
would decrease the allowed rate of environment decoher- 
ence by a similar factor. But dynamical gates are not 
perfect in practice, and the holonomic approach may be 
advantageous if it leads to higher operational precision. 

In Ref. [n| we show that for the BS code our scheme 
can be implemented with at most 3-local Hamiltonians. 
This is optimal for the construction we presented, since 
there are no codes correcting arbitrary single-qubit errors 
with stabilizer or gauge-group elements of weight smaller 
than 2, covering all qubits. An interesting open ques- 
tion is whether it is possible to modify our construction 
so that it uses 2-local interactions, perhaps using recent 
perturbative "gadget" ideas [20| . 

Our abstract scheme proves that the holonomic quan- 
tum computing approach is scalable under a reasonable 



noise model. It is meant as a proof of principle, and 
will undoubtedly require modifications if applied to ac- 
tual physical systems. Given that simple QECCs and 
two-qubit geometric transformations have been realized 
using NMR [2l[ and ion-trap Q techniques, these sys- 
tems seem particularly suitable for hybrid HQC-QEC im- 
plementations. 
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